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ABSTRACT 


The time-averagad power-density distribution below 
the ocean surface due to incident laser radiation is 
examined Бу means of computer simulation of the 
geometrical optics involved with the air/sea interface 
and subsequent ocean penetration by the laser bean. 
The erfects over the entire spectra of incidence 
angles, wind velocities, wind directions, beam spot 


Sizes and depths of penetration are analyzed. 
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I. INTRODUCTION 


a ae m un ee — < _ =F 


The time-averaged power-density distribution below the 
ocean surface due to incident laser radiation is examined 
over the entire spectra of incidence angles, wind velocities 


and wind directions by means of geometrical optics. 


current interest їп using laser systems for 
communications through the air-sea interface and for 
detecting submerged objects necessitates the development of 
prediction rethods for determining the subsurface 
power-distrikution of radiation from a laser source above 
the ocean. A mathematical model F Or the 
optical-communications application done by Karp in Ref. 1 is 
developed in terns of a radiance function related to the 
mutual coherence function With accurate results for 
incidence angles out to 45 degrees. A more general model 
developed Гу Swennen in Refs. 2-4 examines the power-densitjy 
distribution below the ocean surface through a rigorous 
assessment cf the ocean surface geometry allowing for 
theoretically accurate results over the entire incidence 


angle spectrum. 


The computer code developed for this analysis of the 
power-density distribution below the ocean surface is an 


expansion of Swennen's theory[Refs. 2-4]. 





II. DESCRIPTION OF PHYSICAL PARAMETERS 


A. OCEAN SURFACE 


A gecmetrical representation of the ocean surface in 
terms of local slopes developed by Cox and Munk[Ref. 5) is 
the basis for the analysis of the power probability 
distribution kelow the ocean, as it is in Refs. 2-4. A 


review of this representation follows. 


The center of symmetry of the incident laser beam cn tne 
ocean surface facet is the center of a right- handed 
e teslan cocrdinate system, point "О" “in Fig. 1. The 
Z-axis points vertically upward. The y-axis is in the 
horizontal plane, colinear with the projection of the center 
of the incident laser Deam onto that plane and pointing in 


the direction cf the laser source. 


The slopes of the ocean surface facet are defined py 
angular parameters Alpha and Beta. The angle Beta is the 
aide mecetwect the dine of steepest ascent of the facet and 
the x-y plane; Alpha is the angle between the projecticn of 
the line cf steepest ascent onto the x-y plane and the 


y-axis. 


A second cartesian coordinate system at the surface is 
used, as explained in Ref. 2, to simplify the computation of 
pie meslore ss crobability function, This coordinate system is 
designated Ly primes with the y'-axis pointing towards the 


wind source and rotated an angle Chi from the y-axis. The 








z'-axis coincides with the z-axis. The slope parameters are 


also primed with their transformations given by: 


Peta’ = Beta 


Alpha’ = Alpha - Chi 


Belcw the ocean surface a depth Z at the point of 
observation, O', is centered a translational transformation 
of the surface coordinate system. Angles Mu and Nu define 
the laser keam refracted to O! by any surface facet. The 
first angle, Mu, is the angle between the refracted ray na 
the z-axis; the second, Nu, is the angle between the y-axis 


and the vertical plane of the refracted ray. 


B. LASER BEAM 


A cylindrical laser beam of radius a is incident on the 
Ocean surface wave facet at an angle Psi measured from the 
г ахас. The intersection of the beam and the facet is a 
БИЕК Иер саса surface as illustrated in Pig. 2. The 
semi-minor axis of the ellipse is equal to the beam radius 
and the semi-major axis, b, is defined by: 

\ 


b = а * sec Psi 


The analytical development, as in Refs. 2-4, begins at 
the ocean surface facet. The laser source, aining, 
atmospheric propagation, etc. ale extraneous to this 
analysis. Pcwer density values at the surface are not 


quoted because the analysis deals with the ratio of the 








power density at the depth of interest, Eat and the power 


density at the surface, = 5 
5 


ПРАХ TL MUM ANGLE ОР INCIDENCE 


The angle of incidence with respect to the z-axis, Psi, 
of the laser team is analyzed from zero degrees to a maximum 
of 84 degrees. At greater angles the beam will not be able 
to strike all the wave facets within a wavelength cf tne 
wave due to the wave height. There is also the possibility, 
Tee trated in Fig. 3, that a ray incident at Psi > 84 
degrees could fass through the wave, re-enter the atmosphere 
and then re-enter the ocean; or, depending on the local 
Slopes, the ray could undergo a total internal reflection if 
the Critical Angle were exceeded. The mathematical model of 
Swennenf[ Refs. 2-4] and this thesis do not address this 


Situation. 


This lizitation on the maximum angle of incidence, 


Psi foes dcrlved from Pig. 4 and defined ру: 
max 


1 
Psi = tan (L/2H) 


The significant wave height, H, and the average period, T, 
of the wave are given by equations (1) and (2) as defined by 
Pierson[ Ref. 7]. The average wavelength of the wave is 


given by equation (3) as defined by Hillf{ Ref. 8}. 


е2 2 
B = 2.14x10 * W (1) 











TS OST (24) (2) 
2 
L = gf 72 (3) 
where 
W = wind velocity 
g = gravitational constant 


For wind velocities of 1 to 14 n/sec the maximum angle 


of incidence is 84.2 degrees (+0.4 degrees) for all wind 


velocities; therefore, Psi is limited to 84 degrees. 
| max 
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PET SZERORERSIENSETY TERORBILLEITY INTEGRAL 
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A. GENERAL EQUATION 

The power-density probability distribution = ac the 
point o£ cbservation below the ocean surface (0!) due to the 
entire peam is given by the integral of = оуег ба оспе 
contributing surface facets. As a function of the angular 


parameters (Mu,Nu) the integral as defined in Refs. 2-4 is 


given by: 


E = Јаз ама амо (4) 
wv 


where 


F(Mu,Nu) = P. *T*DTF*cos (WI) *sec (WR) 
5 


ШЕ сс u Sp l) costus 


2 
#*P (Z*, 2°) *tan (Beta) *sec (Beta) *J 
х Y 


and 
Р = power density at the surface facet 


T = Fresnel's transmittance coefficient 
pee = daftuse transuittanee function 


WI = angle of incidence 
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WR = angle o£ refraction 


Mu 
0 


first quadrant angle the refracted center 


ray of the beam makes with the positive 
т=ах1< 


Bir = positive angle the projection of the 


refracted center ray of the beam onto the 
xy-plane makes with the y-axis measured 
from the latter to the former 


F(Z!,Z*) = time-average slope distribution 
X y 


Са ср 


The development of equation (4) is described in great 
detail in Refs. 2-4 and will not be repeated here; however, 
for clarity the major components of equation (4) will be 


discussed. 


\ 


BONUUUSSHEL'S TRANSMITTANCE FUNCTION 


Fresnel's transmittance coefficient (T) is defined as 
the intensity ratio of the transmitted to the incident bean 


energy. 


Bi co mated incident laser beam iS "composed of twe 
polarization components, a transverse magnetic or parallel 


Polarlizatıon =, and a transverse electric or normal 


поа 2а ой (T ). Resolving the refracted beam into its 
1 


components (T Т  ) is a complex task due to rotation of the 
1 


plane of incidence about all three coordinate axis over the 
range of integration. In order to deal with this complexity 
the beam is assumed to be rendered diffuse and unpolarized 


below the océan surface due to scattering, propagation 


12 





direction variations and phase changes. This assumption 


allows for the averaging of T and 7 со obtain the value 
L 


Or T. 


ES EET TZ 
R7 


The equations for T апа » used by Swennen[ Refs. 2-š] 


1 
are: 
| . . 2 
E - [sin(WI)*sin(WBR)/sin(WI*WR) ] (5) 
| | и 2 
m = [sin(WI)*sin (WR) /sin (NI+WR) *cos (NI-WR) ] (6) 


Equations (5) and (6) are in disagreement with those 


developed by Ecrn[ äef. 9) and Fowles[Ref. 10], given by: 


T 
ү 


2 
[ 2*cos (WI) *sin (Wa) /sin (WI+NWR) |) (7) 


2 
[ 2*cos (WI) *sin (WR) //sin (WI+WR) *cos (WI-WR) ] (8) 


T 
// 


Equaticns (5) and (6) will give erroneously high values 
IDE r c nc l transmittance function, especially at higher 
angles of incidence. Because of this, equations (7) and (8) 


are used in this modeling. 


The gecmetry of the ray refraction through the ocean 


surface 1S illustrated in Рід. 5. 
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CORRE F RUSE TRANSMITTANCE FUNCTION 


The diffuse transmittance function (DTF) accounts for 
bean attenuation below the ocean surface due to scattering 


and absorbticn. 


The six-ccnstant DTF developed by Duntley in Ref. 11 is 
Simpliried tc a function of two constants; the backward 
scattering coefficient СЕБИ сла The total "absorbtron 
coefficient (АС). This simplification is valid because of 
the randcm scattering particle orientation in the ocean. 
The DIF is defined by: 


DTF = K/[ (AC+BS) *sinh (K*Z) +K*cosh (K*Z) ] 


where 


1/2 
[ AC* (AC* 2BS) ] 


deptn of the point of observation 


IneNcurrent' terminolgy this method of accounting for 


attenuaticn is a zero-angle forward scattering technique. 


Typical values for BS and AC were used fron 
measurements by Tyler[Ref. 12] and Duntley[ Ref. 13]. These 
represent an average between Ccol and warm ocean waters in 


the blue-green spectrum (4800 angstrom). 


BS 0.065 percent 


- 1 
0.088 m 


АС 
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D. SLOPE PRCEABILITY FUNCTION 


rec mezaverage "slope probability function, P(Z',Z'), 
x Y 
is a statistical distribution of the ocean surface slopes 


developed by Cox and Munk[Ref. 5]. The distribution derived 
frcm the surface geometry is Gaussian, which is then altered 
by a Gram-Charlier series in order to account for slope 


Skewness and peakedness caused by the wind. 
Pipe ЕРЕ (Т ЕСа я) 
x Y 


where 
Xi = standardized crosswind slope component 
Eta = standardized upwind slope component 


W = wind velocity 


The development and the equation for the slope 
ПІБОБаАоли т Еу functicn are covered in detail in Refs. 2,4, and 
5% A clean ocean surface is assumed and the linits of 


applicability placed on the functional parameters are 
adhered tc, namely: 


Күл 2.5 
|Еїа| < 2.5 


W < 14 m/sec 


E. JACOEIAN 


The power-density probability integral, equation (4), 


15 








was first developed over the slope components (Z',Z'). It 
x Y 
was then transformed to a function of the angular parameters 


(Alpha,Beta) and then by the Jacobian (J) to a function of 
(mu; Nu). 


> = (far, 


2,2, 


({ F (Alpha, Beta) dAlpha dBeta 
< б 


II 


-1 
( (F ма, ма) 4 амч ани 
им | 


where 


C, 
й 


д (Мо, Ми) И у (А1рћһа, Вета) 
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TI METHODS OF SOLUTION 


ES zu ii u en a A = =< | چ‎ a aE SP = 


A. FAR ZCNE 


Swennen[ Refs. 2-4] observed that when the depth (2) of 
the point of observation is large compared to the bean 
cross-sectional radius at the surface (Z/a>100) tne 
integrand of equation (4) remains essentially constant 
during the ship Se patshs able pare The integral can then be 


approximated by: 


FU 
ii 


-1 
ГР пи, на) жа dMu dNu 
pv 


ie 


= 
Е (Ми, Ми) *J А Мо аби ; with оша & as 


The area aMuáNu in the (Mu,Nu) plane is approximated by 
an ellipse of area AA, whose development is covered in 


detall in Ref. 3. 


2 


This is the Far Zone approximation and equation (4) 


beccmes: 
2 
P = Р, *kpxtan(Bpeta)*sec *сос5(И1) 
s 


d 


-1 
*sec(WR)*P(Z',P')*J *AA 
x у 
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The basic Far Zone computer program developed іп 
FORTRAN language and run on the IBM 360 system solves for 
the ratio of the power density at the point of observation 


to the power density at the surface facet, Ears . The main 
s 


portion cf the program allows for the entry of any 
combination of variables (wind, depth, beam radius, Psi, 


eu Кас; Lackward scattering coefficient, aksorbtion 
coefficient) and also computes the diffuse transmittance 


function, Fresnsl's transmittance function and the slope 
Пао аю Еу function. Subroutines common to all variaple 
entries are used to compute the angles (Alpha, Beta, WI, 
WR), the Jacctian and A. 


It is possible to investigate the power-density 
probability distribution over a wide range of variables with 
a minimai cf storage requirements and computation time due 
to the use of an IBM System/360 Source Library subroutine, 
NLNSYS, in solving for the coupled (Alpha,Beta) and  (Mu,Nu) 


angle pairs from the simultaneous non-linear equations: 


ec e, [cos (Psi) + K*cos (Beta) ]/n (9) 


eu.) cot (Alpha) = sin (Psi) /K 
*sin (Alpha) *sin (Beta) (10) 


where 
n = relative index of refraction (1.33) 


K = K(Psi,Alpha,Beta,n) 
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В. NEAR ZCNE 


When the depth or the point of observation is not large 
compared to the beam cross-sectional radius at the surface 
(Z/a<100) the Far Zone approximation does not hold and the 
double integration of equation (4) must be carried out. 

Swennen{ Ref. 2] used Simpson's method of numerical 
integration to approximate equation (4). The method was 
subject to Singularity points in the integration limits 
which increased the complexity and time required for the 


computation. 


The double integration can be more accurately and 
relatively simply approximated by using the Gauss Quadrature 
system of sclution, Refs. 14 and 15. The Two-Point 
Gauss-Legendre Quadrature method was used. This method, 
valid for polynomials up to degree three (equation (4) is of 
degree two), consists of first transforming the function 
F(Mu,Nu) into a function F(s,t) whose interval is - 1<s<1 and 
IS < BY letting: 


И = [ (Müu —QüMü )*s+Müúu +Mu 1/2 (11) 
u 1 u 1 
а Nu = [ (Nu -Nu )*t+Nu +Nu ]⁄2 (12) 
u 1 u 1 
and 
dMu = [ (Mu -Mu )/2]Jds (13) 
u it 
амо = ( (Nu -Nu )/2 Jdt (14) 
u j 
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The power-density integral 
Hu Yu 
= - Ге (Mu,Nu) dMu dNu 
E 
is transformed into 


14 
P = M - = 
3 [ ( u: au un. oz ]* IN (St) азга (15) 


The Twc-Point Gauss-Legendre Quadrature approximation 


of equation (13) is given by: 


DL 1 
[ (Wu —Mu ) (Nu -Nu )/4 3 Ж. Z u *P(s ,t ) (16) 
u 1 u 1 TON 


P = є 

d 1:0 1:0 1 

where 

W = 1.0 = Gauss weight factors 

2 

= 

s = £(3) = Gauss roots 

1 


The basic Near Zone computer program developed uses 
this Gauss-Legendre Quadrature method in solving the 
power-density integral coupled with the same subroutines as 
the Far Zone progran. The main program computes the 


integration limits (Mu дао) апа (Ми a after determining 
u u 


whether the surface projection of the point of observation 
is inside or outside tne beamts horizcntal cross-secticnal 
area. The equations for the integration limits are well 
defined in Refs. 2-4. 


The transformations described by equations (11)- (15) 


are then performed and the power-density is solved for using 


eguation (16). AS in the Far Zone program, the Near Zone 


20 





program can take any combination of variables as entries. 


2 





Results for oblique incidence angles (0<Psis84 degrees) 


were obtained for vertical plane cuts at чч = 0 degrees 
over the entire range of variables; however, valid data for 
vertical rlanar sections at He * 0 degrees could not be 
obtained due to values of the standardized crosswind and 


upwind slope components being out of the range of 
applicability for the slope probability function. This 
problem is believed to be inherent in the methca of 
Mathematical analysis of the ccean surface used. Results 
for normal incidence (Psi = 0 degrees) over the entire range 


of variables, including Nur planar variations, were 


obtained. 


The validity check carried out on the slope probability 
function, as discussed in Section III.D, clearly identified 
invalid data. Another source of invalid data occurred at 
higher incidence angles (Psi) where the Simultaneous 
solution of ncn-linear equations (9) and (10) resulted in an 
angle of incidence (WI) greater than 90 degrees while the 
Sloe PEOBaLIlIEY function Still indicated valid. This 


output was also easily identified and eliminated. 


Computaticns with the Near Zone solution method used 
more computer time and storage than did the Far Zone 
solution method. For a computer solution using one value 
for each variable entry, the Near Zone used 36 per cent nore 


computer time and 8 per cent more computer storage than did 


D 








an equivalent Par Zone solution. The higher time 
requirement for the Near Zone solution was due mainly to a 
greater compiling time requirement. The time difference 
between the two solution methods decreased as the numter of 
Variables that were incremented during one computer run 


increased. 


A sample computer run with output is presented in 


Appendix В. 
Aw PRESENTATION OF RESULTS 


Results are presented graphically to simplify 
compariscns cf the wide range of variables. Vertical-rlanar 
Sections cf the power density probability distribution below 
the ocean surface cut through the center of the bean 
cross-section at the horizontal ocean. surface at an angle 
Я are used for the majority o£ the results presented. 
This presentation is in rectangular coordinates with 
ncrmalized power density in decibels on the ordinate and the 


angle А, ln degrees on the abscissa. The normalized power 
| 2 
density is equal to 10*log [ (P *Z )/P ].- This 
10 d ds 


: 2 
normalization with respect to > /2 should produce a 
s 


distribution function that is independent of depth in the 


Far Zone regime, according to Swennen[ Ref. 2-4], because the 
multiple slcres seen in a beam cross-section at the surface 
tend to render the radiation diffuse. The power density of 
diffuse radiation at large depths decreases as the square of 


the depth. When the effects of scattering and absorption 
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are included this independence does not occur, as will be 
discussed later; however, the normalization of the power 
density is used in all regimes for continuity and for 
comparisons with Refs. 2-4. 


COSR EFF FECTS CF WIND VELOCITY 


As wind velocity increases, the power density 
distribution spreads out and the maximum power density 
decreases slightly. This effect is presented in Figs. 6-11 
for various angles of beam incidence with respect to the 
vertical (Psi = 0,30,60 and 84 degrees). Figures 10 and 11 
also show the close correlation between results obtained 
frcm the Near Zone and the Far Zone solution methods for Psi 


= 0 and 40 degrees, respectively. 


The spreading effect on the power density distribution 
is due tc larger refraction angles (WR) caused by higher 
wave facet slopes that occur at higher wind velocities. The 


maximum power density decrease is caused by the spreading. 


For an optical communication or detection system an 
increase in wind velocity would mean that the maximum 
intensity cf the beam that could be focused to a point 
peneath the ocean surface would be reduced, but the width of 
receivable signal radiation or the search width would be 


increased. 


C EFFECTS OF INCIDENCE ANGLE 


As the angle of incidence with respect to the vertical 


(Psi) is increased the maximum power density decreases and 
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the distribution shifts away from the vertical. This effect 
can be seen in Figs. 6-11 and more clearly in Figs. 12 and 
З. Figure 12 is a Par Zone solution at a depth ої 50 
meters for Psi = 0,30,60 and 84 degrees while Fig. 13 is a 
Near Zone solution at a depth of 10 meters for the same 
values of Psi. 


The decrease in maximum power density is a result of 
greater reflectance at the surface facet as the angles of 
incidence (WI) increase with increases in Psi. The shift in 
the center of the distributions away from the vertical is 
caused by the higher angles of refraction (WR) that occur as 


WI increases. 


For an cptical communication or detection system an 
increase in Psi would be expected to lower the maximum bean 
intensity below the ocean surface and to shift the point of 


Maximum intensity away from the point of water entry. 


D. EFFECTS CF WIND DIRECTION 


The effects of wind direction are small compared to the 
wind velocity effects. When the wind direction (Chi) is 
perpendicular to the beam (Chi = 90 degrees) there is 
Slightly less spreading of the power density distribution 
ccapared to Chi = Q or 180 degrees. The maximum power 


density is essentially unaffected by variations in Chi. 
Figure 14 shows the effect of variations in Chi for 
normal incidence and Fig. 15 for oblique incidence (Psi = 40 


degrees). 


Wind direction appears to be of little concern to an 


Optical communication or detection system as related to the 
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power density distribution below the ocean surface. 


ES EFFECTS CF SPOT SIZE 


Figures 16 and 17 show the effect cf increasing the 
incident beam spot radius from 0.1 to 0.5 meters for Psi = 
0 and 40 degrees, respectively. This is done with the 
assumption that the energy per unit area remains constant; 
therefore, the total energy must increase. The power 
density spreads out and increases in intensity uniformly 


over the distribution as the spot size is increased. 


Since the equations used in this solution method deal 


with the ratio of a. , these results simply mean that 
s 


more power in a larger beam will produce a more intense and 


more widely distributed pattern below the ocean surface. 


Peer rkCiS OF DEPTH 


The decrease of power density with increases in depth is 
illustrated in Figs. 18 арӣ 19 for Psi equal to 0 and 40 
degrees, respectively. The decrease is linear for the 


ncrmalized power density as can be seen in Fig. 20. 


Scattering and absorption cause this decrease in tne 
power density. As the depth increases the diffuse 
transmittance function (DTF) rapidly becomes an inverse 
function of the sum of the hyperbolic sine of the depth and 
the hyperbclic cosine of the depth, 
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+ 1 
DTF = £[ (sinhZ+cosh2) |) 


which means the: DTF decreases rapidly with increases іп 
depth and approaches zero in the limit. This agrees with 


the thecry cf Preisendorfer[ Ref. 16]. 


Swennen's prediction of a constant normalized power 
density distribution in the Far Zone regime due to the 
diffusing effects of the ocean in Refs. 2-4 holds only when 
the effects of scattering and absorption are ignored. One 
computer run was made with the DTF set equal to one to 


illustrate this effect in Fig. 20. 


Figures 21 and 22 illustrate the decrease in power 
density with depth, non-dimensionalized with respect to spot 
Size, Z/a, again for Psi equal to 0 and 40 degrees, 


respectively. 


An optical communication or detection system will be 
limited greatly by the power density decrease associated 
with depth. 


G. MAXIMUM EOWER DENSITY 


x The maximum obtainable power density available over the 
spectrum of incidence angles (Psi = 0 - 84 degrees) for a 
given set of conditions is presented in Figs. 23 and 24. 
The maximum power density obtainable does not occur for 
normal incidence, but for a small value of Psi, because of 
the fact that as the wind increases the most probable slope 
is not zero but a small angle. The maximum obtainable power 
density then drops continuously with increases in Psi after 


the peak that occurs around Psi = 10 - 15 degrees. 
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Figure 23 illustrates the decrease in the maximum 
obtainable power density for all values of Psi as the wind 
velccity increases. This is due to the spreading effect on 
the power density distribution caused by the wind that was 
described previously. 


Figure 24 illustrates the independence of the maximum 
power density on the wind direction (Chi) for all values of 
Psi. Chi equal to 270 degrees is not plotted but is 
equivalent to the condition at 90 degrees because of the 
Symmetry of the slope probability function about the wind 
direction. 


The angle with respect to the vertical from the point of 
observation Lelow the ocean surface d) at which the 
maximum power density occurs for each value of Psi is 


plotted in Fig. 23% The plot is independent of wind 
velocity and direction. Examining Figs. 16 - 19 it can be 
seen that Fig. 25 is also independent of depth and spot 


size. Mu for maximum power density increases nearly 
linearly f£rcm approximately 1 to 56 degrees as Psi goes from 
0 to 84 degrees. The increase in us 15 due to the power 
density distribution shifting away from the vertical for 
sing valules of Psi as discussed earlier and 


illustrated in Figs. 12 and 13. 


The obvious consideration with an optical communication 
or detection system is the decrease in maximum power density 
that accompanies higher incidence angles and wind 


velocities. The strong independence of the angle ш” at 
which the maximum power density occurs with respect to the 


wind velocity, wind direction, depth and spot size indicates 


28 








that a satellite based laser navigation system for submarine 


usage may be a possibility. 
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The time-averaged power-density distribution kelow the 
ocean surface due to incident laser radiation has been 
examined through computer simulations in both the Near Zone 
and the Far Zcne regimes. The effects of altering incidence 
angles, wind velocities, wind directions, beam spot sizes 


and depths of penetration have been presented. 


The power-density distribution was found to be highly 
dependent on the angle of beam incidence, the wind velocity 
and the attenuation associated with depth of penetration of 
the laser beam. The fact that the location of tne maximum 
power density was found to be fundamentally dependent on 


only the angle of beam incidence indicates a possibility of 


direction finding capabilities by a submerged receiver. 


Improvements in computer systems and analysis techniques 
have allowed an examination of a much broader range of 
variables than were possible  heretofore. Angles of bean 
incidence up to 84 degrees were considered within the range 
of validity fcr restrictions placed on the simulation. 
Above 84 degrees incidence angle the possible interactions 
at the airysea interface make ап analysis of the 
power-density distribution extremely complex, and 
experimentation will be necessary to give indications of the 


feasibility of laser operations in this area. 


This simulation could be extended to an analysis of much 
larger spot sizes that would be expected from a satellite 
based transmitter. Prettyman and Cermak[Ref. 17] suggested 
that the effects of larger spots could be approximated by 
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mE 


assuming each wave facet independent of its neighbors and 
then integrating over the larger spot area, since the small 
spot represents the outer limit of the effects of the ocean 
surface on a laser beam of greater spot size. There would 
have to be a large number of surface facets contained in the 
spct to insure facet independence and accurate 
approximations, hence a strong dependence on the wavelength 


of the wave associated with wind velocity. 


Àn extension of the simulation to pulsed laser radiation 
is also warrented, since this mode of operation may be 
required for air-to-subsurface laser systems. The effects 
of pulsed radiation occur nearly instantaneously as compared 


to the time-averaged analysis performed in this simulation. 


A further refinement of the Simulation would involve a 
more complex modeling of the bean attenuation due to 
scattering and absorption in the ocean through the use of a 
diffusicn or a multiple forward scattering method, both 
described by Bucherf{Ref. 18] and Preisendorferf[ Ref. 16]. 
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APPENDIX А 


FIGURES 
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APPENDIX B 


COMPUTER PROGRAM USAGE 


TO use the computer prograns listed following this 
Appendix, first compute the ratio of the depth of the point 
of observatıon below the ocean surface to the beam radius at 
the surface facet, Z/a. Use the Near Zone program for 
2/а<100 апа the Far Zone program for Z/a>100. The Near Zone 
program may be used for any value of Z/a, but a slight 


Computation time and computer storage increase will result. 


The only other basic requirement is the selection of 
eight input variables. The variables are listed below along 


With the line from the respective program in which they 


appear. 

NEAR FAR 
VARIABLE ZONE ZONE 
Depth, Z (m) 470 440 
Beam radius, RAD (m) 480 450 
Angle iur NUORUN (deg) _ 490 480 
Backward scattering coefficient, BS 560 530 
Absorpticn coefficient, AC (m ) 570 580 
Beam incidence angle, PSIRUN (deg) 630 590 
Wind direction, CHI (rad) 680 640 
Wind velccity, W (m/sec) 710 670 


Z, RAD, NUORUN, BS and AC are direct program inputs 
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made on the indicated lines. PSI, CHI and W are selected 
from DATA statements 430 and 440 in the Near Zone progran 
and 420 and 430 in the Far Zone progran. The variable 
locations in DATA statements were used for multiple 
Variations of the parameters in a slightly different form of 
the program than presented which is easily set up with DO 
LOOPS. The entry variables may also be varied in this same 


manner. 
Both programs must be combined with the common package 
of subroutines that are listed following the basic program 


listings. 


The cutput from both listed programs appear on the next 
two pages. 
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